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Abstract — Extending network lifetime of battery-operated de- 
vices is a key design issue that allows uninterrupted information 
exchange among distributive nodes in wireless sensor networks. 
(\ Collaborative beamforming (CB) and cooperative transmission (CT) 
C™j have recently emerged as new communication techniques that 
enable and leverage effective resource sharing among collabora- 
iive/cooperative nodes. In this paper, we seek to maximize the 
lifetime of sensor networks by using the new idea that closely located 
n odes can use CB/CT to reduce the load or even avoid packet 
forwarding requests to nodes that have critical battery life. First, we 
study the effectiveness of CB/CT to improve the signal strength at a 
ly-^ faraway destination using energy in nearby nodes. Then, a 2D disk 
case is analyzed to assess the resulting performance improvement. 
For general networks, if information-generation rates are fixed, 
r— ,the new routing problem is formulated as a linear programming 
p—H problem; otherwise, the cost for routing is dynamically adjusted 
KH according to the amount of energy remaining and the effectiveness 
rj^ of CB/CT. From the analysis and simulation results, it is seen that 
(j the proposed schemes can improve the lifetime by about 90% in 
'the 2D disk network and by about 10% in the general networks, 
compared to existing schemes. 
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I. Introduction 

In wireless sensor networks [1], extending the lifetime of 
battery-operated devices is considered a key design issue that 
increases the capability of uninterrupted information exchange 
and alleviates the burden of replenishing batteries. In [2], a data 
routing algorithm has been proposed with an aim to maximize 
the minimum lifetime over all nodes in wireless sensor networks. 
A survey of energy constraints for sensor networks has been 
studied in [3]. In [4], the network lifetime has been maximized 
by employing the accumulative broadcast strategy. The work in 
[5] has considered provisioning additional energy in the existing 
nodes and deploying relays to extend the lifetime. 

Recently, collaborative beamforming (CB) [6] and cooperative 
transmission (CT) [7] [8] have been proposed communication 
techniques that fully utilize spatial diversity and multiuser di- 
versity. While most existing work in this area concentrates on 
improving the performances at the physical layer, CB and CT 
also have impact on the design of higher layer protocols. In 
this paper, we investigate new routing protocols to improve the 
lifetime of wireless sensor networks using these two techniques. 

First, we study the fact that CB/CT can effectively increase the 
signal strength at a destination node, which in turn can increase 
the transmission range. We obtain a closed-form analysis of the 
effectiveness of CT similar to that given for CB in [6]. Then, 
we formulate the problem as a maximization of the network 
lifetime, defined until the time of the first node failure. The new 
idea is that closely located nodes can use CB/CT to reduce the 
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loads or even avoid packet forwarding requests to nodes with 
critical battery lives. From the analysis of a 2D disk case using 
CB/CT, we investigate how battery-depleting nodes close to the 
sink can be bypassed. Then we propose algorithms for a general 
network situation. If the information-generation rates are fixed, 
we can formulate the problem as a linear programming problem. 
Otherwise, we propose a heuristic algorithm to dynamically 
update costs in the routing table according to the remaining 
energy and effectiveness of collaboration. From the analysis 
and simulation results, the proposed new routing schemes can 
improve the lifetime by about 90% in the 2D disk network 
compared to the pure packet forwarding scheme, and by about 
10% in general networks, compared to the schemes in [2]. 

This paper is organized as follows: In Section [II] the system 
model is given, and the abilities of CB/CT to enhance the des- 
tination signal strength are studied. In Section [III] we formulate 
the lifetime maximization problem, analyze a 2D disk case, and 
propose algorithms for general network situations. Simulation 
results are given in Section [IV] and conclusions are drawn in 
Section [V] 

II. System Model and Effectiveness of CB/CT 

We assume that a group of sensors is uniformly distributed 
with a density of p. Each node is equipped with a single ideal 
isotropic antenna. There is no power control for each node, 
i.e., the node transmits with power either P or 0. There is no 
reflection or scattering of the signal. Thus, there is no multipath 
fading or shadowing. The nodes are sufficiently separated that 
any mutual coupling effects among the antennas of different 
nodes are negligible. 

For traditional direct transmission, a node tries to reach 
another node at a distance of A. The signal to noise ratio (SNR) 
is given by 

PC Q A~ a \h\ 2 



r 



(i) 



where Cq is a constant that incorporates effects such as antenna 
gains, a is the propagation loss factor, h is the channel gain, and 
a 2 is the thermal noise level. We define the energy cost of such 
a transmission for each packet to be one unit. 

In Figure Q] we show the system model with CB/CT. In 
traditional sensor networks, the only choice a node has is to 
forward its information toward the sink. This will deplete the 
energy of the nodes near the sink, since they have to transmit 
many other nodes' packets. To overcome this problem, in this 
paper, we propose another choice for a node consisting of 
forming CB/CT with the nearby nodes so as to transmit further 
towards the sink. By doing this, we can balance the energy usage 
of the nodes having different locations and different remaining 
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Fig. 1: System Model 

energy. In the rest of this section, we study how effectively CB 
and CT can improve the link quality. 

A. Effectiveness of Collaborative Beamforming 

Suppose there are a total of N users for collaborative beam- 
forming within a disk of radius R. We have 

N = [pirR 2 \ ■ (2) 
Each node has polar coordinate (rk,ipk) to the disk center. 
The distance from the center to the beamforming destination 
is A. The Euclidean distance between the k th node and the 
beamforming destination can be written as: 



(3) 



where (f> is azimuth direction and is assumed to be a constant. By 
using loop control or the Global Positioning System, the initial 
phase of node k can be set to 

</>fc = -y«), (4) 

where A is the wavelength of the radio frequency carrier. 
Define z = [zi . . . zn] t with 

z k = T7 sinOfc - 0/2). 
H 

The array factor of CB can be written as 

N 



(5) 



F(<j>\z) = — > e 



N ^ 
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The far-field beam pattern can be defined as: 

P(0|z) = |F(0|z)| 2 



(6) 



(7) 
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where 



, „ AnR sin % 

o-W) = — 



(8) 



Define the directional gain D^b as the ratio of radiated 
concentrated energy in the desired direction over that of a single 
isotropic antenna. From Theorem 1 in [6], for large ^ and N, 
the following lower bound for far-field beamforming is tightly 
held: 

D CB 1 

(9) 
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where p « 0.09332. 

Considering this directional gain, we can improve the direct 
transmission by a factor of . Notice that this transmission 
distance gain for one transmission is obtained at the expense of 
consuming a total power of N units from the nearby nodes. 



B. Effectiveness of Cooperative Transmission 

Similar to the CB case, we assume N users are uniformly 
distributed over a radius of R. The probability density function 
of the users' radial coordinate r is given by 



q(r) 



2r 

Wr 2 ' 



0<r<R, 



(10) 



and the users' angular coordinate ip is uniformly distributed 
between [0, 2ir). 

Suppose at the first stage, node 1 transmits to the next hop or 
sink. Then in the following stages, node 2 to node N relay the 
node l's information if they decode it correctly. The received 
signals at node 2 to node k at stage 1 can be expressed as: 



Zk 



Pr^hlx + nl, k = 2, 



.N, 



(11) 



and the received signals at the destination in the following stages 
are 



Uk = J Pd k a hkX + n k . 



(12) 



Here P is the transmitted power, h r k and hk are the channel 
gains of source-relay and relay-destination, which are modeled as 
independent zero mean circularly symmetric complex Gaussian 
random variables with unit variance, x is the transmitted data 
having unit power, and nk and n r k are independent thermal noises 
with noise variance a 2 . 

Theorem 1: Define D^J to be the energy enhancement at the 
destination node due to CB. Under the far-field condition and 
the assumption that channel links between source and relays are 
sufficiently good, we have the following approximation: 



D 



CT 1 + (JV-1) ^(2 -L;a±g;^ 
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(13) 



where L is the frame length and 2F1 is the Hypergeometric 
function 

^ (a) n (b) n z n 



2 Fi(a, 6; c; z) 



n ( C )« n ' ' 



(14) 



where (a)„ = a(a+l) • • • (a+n— 1) is the Pochhammer symbol. 

Proof: The SNR received by the k th user at stage one can 
be written as 



Tk = 



(15) 



where |/i£| 2 is the magnitude square of the channel fade and 
follows an exponential distribution with unit mean. 

Without loss of generality, we suppose that BPSK modulation 
is used and Cq — 1. The probability of successful transmission 
of the packet with length L is given by: 

For fixed (rk,ipk), the average energy that arrives at the 
destination can be written as: 



N 

= J2 Pd k aP r 

fe=l 



(17) 
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Fig. 2: Approximation of CB Effectiveness 

Since for node 1, r\ — 0. We can write the average energy gain 
in the following generalized form: 



N 

D CT = V 

av 

k=l 



R r 2-K 



A a 
2^ 



d^ a P r k (r k )q(r k )dr k d^ k . 



(18) 



Since each user is independent of the others, we omit the 
notation k and can rewrite (f]~8T ) as: 



D 



CT 



l + (N-l) (19) 
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A~ a R 2 



(A 2 + r 2 - 2rAcosip)~2 P r (r)drdip. 



With the far field assumption, we have 



(A 2 +r 2 - 2rAcosip)-%di> « A~ a . 



The average energy gain is approximated by 
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l + (N-l)j^J rP r (r)dr. 



(20) 



(21) 



With the assumption of sufficiently good channels between 
sources and relays, we have the following approximation of ( fTol i: 



P r (r) w (1 



4P 



-) 1 



(22) 



Since 



Hi - T5-) dr = o R 2^1 -,-L; — 
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(23) 

we can obtain ( fT3l . ■ 
In Figure [2] we compare the numerical and analytical results 
of D av for different radii R. Here A = 1000m, P = lOdbm, 
cr 2 = -70dbm, a = 4, L = 100 and TV = 10. We can see that 
the numerical result fits the analysis very well, which suggests 
that the approximation in ( fT3l is a good one. 

III. CB/CT Lifetime Maximization 

In this section, we first define the lifetime of sensor networks 
and formulate the corresponding optimization problem. Then, by 
using a 2D disk case, we demonstrate analytically the effective- 
ness of lifetime saving using CB/CT. Finally, two algorithms are 
proposed for general network configurations. 
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Fig. 3: CB/CT Routing Model 

A. Problem Formulation 

In Figure [3] we show the routing model with CB/CT A 
wireless sensor network can be modeled as a directed graph 
G(M, A), where M is the set of all nodes and A is the set of 
all links £ M. Here the link can be either a direct 

transmission link or a link with CB/CT. Let Si be the set of 
nodes that the i th node can reach by direct transmission. Denote 
by C|™ the set of nodes that node i needs to apply CB/CT with in 
order to reach node m. In the example in Figure [3] C|" = {1} 
and Cf = {}. A set of origin nodes O where information is 
generated at node i with rate Qi can be written as: 



O = {i\Qi > 0, i £ M}. 
A set of destination nodes is defined as D where 

D = {i\Qi < 0, i £ M}. 



(24) 



(25) 



Define q = {% } to represent the routing and the transmission 
rate. There are many types of definitions for lifetime of sensor 
networks. The most common ones are the first node failure, the 
average lifetime, and a lifetime. In this paper, we use the lifetime 
until first node failure as an example. Other types of lifetime can 
be examined in a similar way. Suppose node i has remaining 
energy of Ef. The lifetime for each node can be written as: 

Ti (q) = ^ , (26) 

where the first term in the denominator is for direct transmission 
and the second term in the denominator is for CB/CT. Notice 
that CJ 1 is not a function of q. We formulate the problem as 

max minT; (27) 
q 



s.t. 



qij > 0,Vi, j 



where the second constraint is for flow conservation. 

B. 2D Disk Case Analysis 

In this subsection, we study a 2D disk case network. Users 
with the same remaining energy are uniformly located within a 
circle of radius Bq- One sink is located at the center location 
(0,0). Each node has a unit amount of information to transmit. 
Here we assume the user density is large enough, so that each 
node can find enough nearby nodes to form CB/CT to reach the 
faraway node. 
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TABLE I: Lifetime Saving vs. Disk Size 



Fig. 4: Analytical Results for the 2D Disk Case 

For traditional packet forwarding without CB/CT, the number 
of packets needing transmission for each node at the distance B 
to the sink is given by 

N pf (B)= J2 



1 



nAr 



D 



(28) 



where Aq is the maximal distance over which a minimal link 
quality 70 can be maintained, i.e. T(Aq) — 70. 

If all nodes use their neighbor nodes to communicate with the 
sink directly, we call this scheme pure CB/CT. To achieve the 
range of B, we need Nqb/ct{B) for CB/CT, i.e., 



D„ 



1 N C b/ct{B ) 
pn 



(max(|-,i; 



For collaborative beamforming, we can calculate 



Ncb(B) > ~ ( <uC2 



(29) 



(30) 



where cq = (max(j-, 1)) Q and c\ = n\y/]m. For cooperative 
transmission, numerical results need to be used to obtain the 
inverse of 2-F1 in Theorem 1. Notice that if the node density is 
large enough, then D av /N — > 1. 

In Figure@] we show the average transmission per node vs. the 
disk size B^. We can see that for traditional packet forwarding, 
the node closest to the sink has the most transmissions per node, 
i.e., it has the lowest lifetime if the initial energy is the same for 
all nodes. On the other hand, for the pure CB/CT scheme, more 
nodes need to transmit to reach the sink directly when Bq is 
larger. The transmission is less efficient than packet forwarding, 
since the propagation loss factor a is larger than 1. The above 
facts motivate the joint optimization case where nodes transmit 
packets with different probabilities over traditional packet for- 
warding and CB/CT. 

For traditional packet forwarding, nodes near the sink have 
lower lifetimes. If the faraway nodes can form CB/CT to transmit 
directly to the sink and bypass these life depleting nodes, the 
overall network lifetime can be improved. Notice that in this 
special case, if the faraway nodes form CB/CT to transmit to 
nodes other than the sink, the lifetime will not be improved. 
For each node with distance B to the sink, and supposing the 
probability of using CB/CT is P r (B), we have 
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P r (B + jAo)), 



(31) 



where the first term on the right-hand side (RHS) is the necessary 
energy for transmitting one packet, and the second term is the 
number of packets for transmission. The goal is to adjust P r (B) 
such that the lifetime is maximized, i.e., 



min max N. 



joint 



(B). 



(32) 



l>P r (S)>0 

Notice that N CB / CT > 1, and in ( fJTb the second term on the 
RHS depends on the probabilities of CB/CT being larger than 

B. So we can develop an efficient bisection search method to 
calculate (l32l . We define a temperature n that is assumed to 
be equal or greater than Nj i nt (B),\/B. We can first calculate 
Nj i n t(B) from the boundary of the network where the second 
term on the RHS of OTb is one. Then we can derive all 
Nj i n t(B) by reducing B. If k is too large, most information 
is transmitted by CB/CT, and the nodes faraway from the sink 
waste too much power for CB/CT; on the other hand, if k is 
too small, the nodes close to the sink must forward too many 
packets. A bisection search method can find the optimal values 
of k and N joint (B). 

In Figure 2J we show the joint optimization case where the 
node density is sufficiently large. We can see that to reduce the 
packet forwarding burdens of the nodes near the sink, the faraway 
nodes form CB/CT to transmit to the sink directly. This will 
increase the number of transmissions per node for them, but 
reduce the transmissions per node for the nodes near the sink. 
In Table U we show the maximal Nj i n t{B) and the lifetime 
saving over the traditional packet forwarding. We can see that 
the power saving is around 90%. 

C. General Case Algorithms 

In this section, we first consider the case in which the 
information generation rates are fixed for all sensors, and develop 
a linear programming method to calculate the routing table. Here 
to simplify the calculation of set C™, we assume its size equals 
one. Obviously, this is suboptimal for ( T27l i. Then we select 
the nearest neighbor for CB/CT. C™ = 1, if node j's nearest 
neighbor can help node j to reach node m. Define q^ = Tq^. 
The problem can be written as a linear programming problem: 

max T (33) 

Qij > 0,Vi,j; 

where the second constraint is the energy constraint and the third 
constraint is for flow conservation. 

Next, if the information rate is random, each sensor dynam- 
ically updates its cost according to its remaining energy and 
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Fig. 5: Snapshot of CB/CT Routing 



Fig. 6: Lifetime Improvement Comparison 



with consideration of CB/CT. Some heuristic algorithms can be 
proposed to update the link cost dynamically. Here the initial 
energy is Ei. Define the current remaining energy as E_ i . We 
define the cost for node i to communicate with node j as 



cost, 



01 



iec J 4 



(34) 



where f3i and 02 are positive constants. Their values determine 
how the packets are allocated between the energy sufficient and 
energy depleting nodes, and between the direct transmission and 
CB/CT. Notice that ( f34b can be viewed as an inverse barrier 
function for E_ t > 0. 

IV. Simulation Results 

We assume nodes and one sink are randomly located within 
a square of size LxL Each node has power of lOdbm and 
the noise level is -70dbm. The propagation loss factor is 4. The 
minimal link SNR is lOdB. The initial energy of all users is 
assumed to be unit and information rates for all users are 1 . 

In Figure |5J we show a snapshot of a network of 5 sensor 
nodes and a sink with L = 50m. Here node 1 is the sink. The 
solid lines are the links for the direct transmission, and the dotted 
line from node 6 to the sink is the CB/CT link with the help 
of node 5. For traditional direct packet forwarding scheme, the 
best flow is 
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with the resulting energy consumed for all nodes given by 
[0,1.0,0.2,0.2,0.3,0.3]. Because node 2 is the only node that 
can communicate with the sink, the best lifetime of this routing 
is 0.2 before node 2 runs out of energy. 
With CB/CT, the best flow is 
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with the energy consumed for all nodes given by 
[0,1.000,0.333,0.333,1.0000,0.782]. Here some flow can 
be sent to the sink via node 6. Because of CB/CT, node 5 has 
to consume its power. The lifetime becomes 0.333 which is 
67% improvement over direct packet forwarding. 

In Figure [6] we compare the performance of three algorithms, 
the shortest path, the algorithm in [2], and the proposed CB/CT 
algorithm. Here L = 100m. As the number of users increases, 
the performance of the shortest path algorithm decreases. This 
is because more users will need packet forwarding by the nodes 
near the sink. Consequently, they die more quickly. Compared 
with the algorithm in [2], the proposed schemes have about 10% 
performance improvement. This is because of the alternative 
routes to the sink that can be found by CB/CT. 

V. Conclusions 

In this paper, we have studied the impact of CB/CT on 
the design of higher level routing protocols. Specifically, using 
CB/CT, we have proposed a new idea based on bypassing energy 
depleting nodes that might otherwise forward packets to the sink, 
in order to improve the lifetime of wireless sensor networks. 
From the analytical and simulation results, we have seen that 
the proposed protocols can increase lifetime by about 90% in 
a 2D disk case and about 10% in general network situations, 
compared with existing techniques. 
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